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In this work we study the process of depinning of a quantum of circulation trapped inside a disk 
by an apphed two dimensional superflow. We use the Gross-Pitaevskii model to describe the neutral 
superfluid. The collective coordinate dynamics is derived directly from the condensate equation of 
motion, the nonlinear Schrodinger equation, and it is used to obtain an expression for the critical 
velocity as a function of the defect radius. This expression is compared with a numerical result 
obtained from the time independent nonlinear Schrodinger equation. Below the critical velocity, we 
obtain the dependence of the semiclassical nucleation rate with the flow velocity at infinity. Above 
the critical velocity, the classical vortex depinning is illustrated with a numerical simulation of the 

S ' time dependent nonlinear Schrodinger equation. 
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I. INTRODUCTION 



A superflow passing through a constriction or an obstacle of some sort can produce vorticity in a process known as 
inhomogeneous nucleation [Q|j . The vortex is created at the boundary between the superfluid and the container walls 
■ O ■ and two flow regimes are identified. When the flow velocity, far away from the walls, is smaller than a critical velocity 
C ' an effective barrier is present that prevents vorticity to appear in the flow, this is the laminar flow regime. Above the 
critical velocity, a vortexless state can dynamically produce vorticity. In the laminar flow regime, and in the absence 
of thermal excitations, it is believed that vorticity can still be produced via Macroscopic Quantum Tunneling (MQT). 
Measurements of the vortexproduction rates for a sphere moving in a superfluid and for an oscillating superflow 
through a microscopic orifice |Q,|| have shown evidence of MQT of vorticity at sufficiently low temperatures. Numerical 
studies of the dynamical nucleation process, above the critical velocity, have relied on the nonlinear Schrodinger model 
for the superfluid dynamics 1^,0. Analytical studies of the tunneling process have used the collective coordinate 
^~^ ' ' approach to map what is really a fleld theory tunneling problem into a standard Quantum Mechanical problem for 
[ — , the vortex coordinates |^. Besides the bosonic superfluid, similar ideas have been applied to vortex states in short 
^> ' coherence length superconductors ^-p^ . 

'"^ Different from Ref. [@| , that studied the classical nucleation of a vortcx/antivortex pair in the proximity of a circular 

^^ ■ obstacle in a two dimensional superflow, in this work we will investigate the depinning of a quantum of circulation 
trapped inside the obstacle. 

In the following sections we define the Gross-Pitaevskii model and extract the dynamics of the collective coordinate 
from the condensate equations of motion. We obtain an expression for the critical velocity and compare it with a 
i-j^ . numerical result. The process of dynamical nucleation, above the critical velocity, is illustrated with a numerical 
(~| ' integration of the nonlinear Schrodinger equation. The Macroscopic Quantum Tunneling, below the critical velocity, 
Q . is studied within the collective coordinate approach and the semiclassical nucleation rate is obtained. 
O 

> 

C^ ' The problem we consider is that of a two dimensional superflow past an impenetrable disk perpendicular to the 

z axis. We will be interested in the process of removal of a preexistent quantum of vorticity inside the disk by the 
applied superflow (along the y direction). In three dimensions this would correspond to a straight vortex depinning 
from a columnar defect due to an applied superflow perpendicular to the column axis. In Ref. [[l2[ this problem is 
studied considering that the vortex depins in the form of an arc attached to the column. 

We start from the Gross-Pitaevskii model for the dynamics of the condensate wave function, i.e. the nonlinear 
Schrodinger equation (NLSE) |l|,|l|: 

Bill fi^ 

in^ = -—V^ + g{\^\^-po)^. (1) 



II. THE MODEL 



Here m is the boson mass, g is the parameter in the hard core potential between the bosons and po is the equihb- 
rium superfluid number density. The obstacle and any applied flow impose appropriate boundary conditions on the 
condensate wave function if). 



It will prove convenient to express length in units of the coherence length C — \/^ /ffPo^^i the condensate in units 
of ^/p^ and velocities in units of the speed of sound c = fi/niS,. 
We take as boundary conditions on ?/;: 

^(r,t)^e™~^e*^ r ^ oo, (2) 

V'(r,i)=0, r<a. (3) 

The first and second exponentials in y say that at infinity we have the external flow along the y direction and the 
circulation due to a vortex somewhere near, or inside the disk. We assigned a radius a to the disk and placed the 
origin at its center. 

The determination of the critical velocity for vortex depinning can be made by looking for steady state solutions 
of the NLSE equation above. For Voo below a certain critical value, dependent only on the disk radius, one is able 
to find stationary solutions representing a laminar flow. Above the critical velocity no stationary state exists, a state 
with a circulation around the disk can evolve to produce a vortex near the edge of the disk. The NLSE in Id was 



found to display the same type of behavior |15 lq|, phase slips are dynamically produced behind an inhomogeneity if 



a large enough flow velocity is applied at infinity. We postpone the discussion of the numerical calculation to section 



IV. 



III. THE COLLECTIVE COORDINATE APPROACH 

Now we will follow Ref. B and study this problem using as a collective coordinate the 2d location of the vortex in 
the flow region. We will ignore the role of the superfluid compressibility. 

When we write "0 = ^/p e*"^, the NLSE yields a continuity equation and an Euler equation for the superfluid density 
p and superfluid velocity v = V0. In the incompressible limit these equations become: 

V • V = 0, (4) 

|^ + (v.V)v = 0. (5) 

If a set of vortices of charge rii = ±1 is present at positions R,;, we also have: 

V X V = ^ 2Tm,5{r - R,). (6) 

i 

The boundary condition |3| becomes 

V • f = 0, r = a. (7) 

The problem of determining the velocity field from a given distribution of vortices and a fixed boundary condition 
at infinity is contained in equations and l6l They can be cast in the form of a 2d electrostatic problem if we define 
the potential x'- 

V = Vx X z. (8) 

Eqs. and H reduce to a Poisson's equation: 

V'x = -5I2^"«'^(r-R^). (9) 

i 

As for the boundary condition, Eq. translates to x ^^ Xo ^ WooS^ at infinity and the disk becomes an equipotcntial 
of X due to the impenetrability condition [7|. 



The vortex dynamics is contained in equation 0. As is well known, inside an ideal and incompressible fluid the 
vortex follows the flow velocity at its center (excluding the singular contribution of the vortex itself to v) . This being 
simply Kelvin's theorem p7|| : 

^=v(R,i). (10) 

In the Appendix we show how to go from Eq. || to Eq. \1(] . 

We now derive the dynamics of a vortex of charge n at position R outside the disk, with an applied flow at inflnity, 
Voo — Vooilj by first writing the general solution of Eq. ^ as the sum of three terms: 

Xc(r) = -Wooa:(l-aVr2), (11) 

Go(r,R)==-nln|r-R|, (12) 

G/(r, R) = +nln |r - {a^/R^)Ii\ - n ln(r/a). (13) 

The total flow, which obeys the topological constraint § w ■ d\ = n for a path enclosing the disk at infinity, is given 
by: 

v(r, t) = Vr [xe(r) + Go(r, R) + G/(r, R)] x z. (14) 

To find the flow velocity at the vortex position and the vortex equation of motion we use the symmetry properties 
of the Green's functions and obtain (excluding the singular contribution of Gq): 






Xe(R) + ^Gz(R,R) 



X z. (15) 



One can view this equation as Newton's equation for a massless object moving in a constant _B-field along the z 
direction plus an electrostatic potential set by Xc and G/ (the applied flow plus the vortex images). Alternatively, the 
X and Y coordinates of the vortex can be thought of as conjugate variables with Xc + \Gi playing the role of the 
Hamiltonian. 

The absence of a mass term is due to the incompressibility assumption, the effect of the compressibility is to 
introduce retardation effects in the response of the vortex to its external surroundings, which in turn gives rise to 
a frequency dependent mass term in the vortex equation of motion |l8[ . The issue of the vortex mass in a neutral 
superfluid was, until recently, the subject of many papers p9| ]. 

Fig. |l| shows the isolines of the Hamiltonian function, which correspond to the trajectories followed by a positive 
vortex in the incompressible limit. There one sees a saddle point to the right of the disk separating closed trajectories 
near the disk from open ones, of the same energy, further from it. 

As it stands. Fig. H describes the dynamics of a vortex that is already outside the disk, but we want to describe 
the process of vortex depinning whereby the quantum of circulation, initially inside the disk, is pushed outwards by 
the flow in the form of a vortex. The problem we face is characteristic of the collective coordinate approach since we 
are using as variable the very object whose creation we want to describe. One way out of this dilemma is to recognize 
that the incompressibility hypothesis should hold only at a distance 5{'^ ^) away from the obstacle. The nucleation 
process in fact happens, due to compressibility effects, within a coherence length from the disk, see Ref. 0] and Fig. 
below. 

If the parameter v^o is such that the saddle point is inside this region dominated by the compressibility, one has 
the possibility of the nucleated vortex to appear in one of the open orbits to the right of the disk and be dragged 
along by the flow, as one observes in the simulation of the NLSE. Correspondingly, when the saddle point is outside 
the compressibility region any nucleated vortex would remain close to the disk, in one of the closed orbits. We would 
like to identify this situation with what happens when the velocity is below the critical velocity and the quantum of 
circulation remains trapped inside the disk. 

We therefore flnd an estimate for the critical velocity from: 



d 
dX 



Xc(R) + ^G/(R,R) 



= 0, (16) 

{X,Y} = {a+Sfi) 



and view S as an adjustable parameter of order ^. 



We obtain: 



[a + .5)4 - a4 ■ 



(17) 



An alternative expression for the critical velocity is that used in Ref. [gO[. One writes the NLSE as two equations 
for p and 0. One of the equations is a continuity equation, the other is: 

In the steady state, d<j)/dt must be a constant. The asymptotic region fixes this value and one finds: 

If one ignores the Laplacian in this equation one obtains a relation between v and p. Above a certain value of v 
one has d[p{v)v]/dv < 0, this defines Wmax — \/(w^ + 2)/3 as the maximum possible velocity inside the fluid. If at 
some point v > Wmax a vortex will be nucleated so as to reduce the flow velocity at that point. We can estimate the 
highest velocity in the flow using the incompressible limit. If a quantum of circulation is present inside the disk, v is 
maximum at the equator, and it is equal to 2t;oo + 1/a. This gives the critical velocity as the maximum value of v^o 
that ensures v < fmax everywhere in the flow: 



3 

11 



2 1 22 



(20) 



In the limit a ^ oo one finds the critical velocity for vortex/ antivortex pair nucleation, (wc)pair = ■\/2/ll, quoted 
in Ref. §. 

In Fig. ^ we show the dependence of Vc on the disk radius from Eqs. |l^ and ^. In the first case we used d — 0.5. 

IV. NUMERICAL SIMULATION 

To determine numerically the critical velocity for depinning we looked for the values of Uoo which allowed a steady 
state solution of Eq. |l|. This solution, with the boundary condition of Eq. at infinity, has e~'^°°*/^ as time 
dependence. The spatial dependence is then found as a solution of 

^V = -^vV + (I^P-i)V', (21) 

with boundary conditions given by Eqs. and S. 

We wrote this equation using finite differences on a two dimensional grid and reduced the problem to a set of N 
complex, nonlinear algebraic equations, for N complex variables, the values of ip at the grid nodes. This was solved, 
starting from a trial state, using Newton's method ^^. The number of nodes varied depending on the disk radius so 
that we always had a minimum distance of 10^ between the disk boundary and the sides of the working area. The 
grid spacing used was 0.5^, but it was reduced to 0.25^ for disks of smaller radii to ensure a reasonable representation 
of the disk as a region of the square grid. The disk radii ranged from 0.5^ to 10^. 

After a solution with a given value of tioo was found, a new Uqo, increased by 0.1, was used in the boundary condition 
and in equation p^. The code was run again taking the available solution as trial state. When Newton's method failed 
to converge we took the last value of Uoo as a lower limit for the critical velocity. 

The critical velocities so obtained are shown in Fig. ||. There one sees an increase of Vc with the disk radius, as 
both Eqs. O and Ofl predict, and an apparent saturation at the value 0.5c, higher than the prediction of y'2/11 c of 
Eq. |2y. The value of S that must be used in Eq. ^to obtain this asymptotic value of v^ is 0.5^. 

Above the critical velocity, if one uses as initial state a condensate with a quantum of circulation inside the disk, 
one can observe the creation of a vortex from the dynamics of the NLSE, Eq. |l|. This equation was solved numerically 
using finite differences in a grid similar to the one described above. The condensate at time zero was evolved using 
a split-step algorithm |2^] where the Laplacian part and the nonlinear part of the NLSE were used alternately. The 
time step was 0.1^/c. The condensate amplitude at four different times is shown in Fig. |[ There one observes the 
vortex formation in the region close to the disk equator and its subsequent movement downstream. After the vortex 
depinning, a vortex/antivortex pair is seen nucleating as in Ref. H. This happens because the value of Woo = 0.5c 
used is also above the critical velocity for pair nucleation for that particular radius, a — 1^. 



V. MACROSCOPIC QUANTUM TUNNELING 



Below the critical velocity an effective barrier prevents the circulation around the disk to become a vortex outside 
the disk. At zero temperature this barrier can only be overcome via Macroscopic Quantum Tunneling. In a thorough 
semiclassical treatment we would have to find the bounce in imaginary time that connects the condensate correspond- 
ing to the stationary, laminar flow, solution of the NLSE at t = — cx) to a condensate containing a nucleated vortex 
near the disk at r = 0, and back to the stationary state at r = +00. This bounce would be a nontrivial solution of 
the imaginary time NLSE [ll^ : 



OT 2 



5^ 
9^ 



1 



VV + (-0V-l)^, V'(r, +00) ^^tZ-oM. 



(22) 
(23) 



Here, "00(1") is the stationary condensate that represents the laminar flow around the disk with the circulation 
trapped inside. The bounce solution would show, as imaginary time progresses, a vortex leaving the disk towards the 
flow and then returning back. The vortex position at t = would indicate the most favorable position inside the flow 
for the vortex to emerge. 

Instead of retaining all degrees of freedom of the condensate in the calculation of the bounce one can use the 
collective coordinate approach of section [II and focus on the dynamics of the vortex position. In this approach the 
bounce, or tunneling path, connects a closed orbit in the XY plane to its open counterpart having the same value of 
E ^ Xc + \Gi. This is done by allowing the momentum Y to become imaginary, see Fig. |[ The logarithm of the 
WKB tunneling rate is then proportional to the pq part of the action that governs the vortex dynamics, Eq. ^. This 
action in real time is: 



The semiclassical tunneling rate is then: 



EwKB oc exp 



yX-Xc(R)-^G/(R,R) 



-2TTpoe f\Y{X,E)\dX 



dt. 



with Y{X,E) being the trajectory in the XY plane that is a solution (with imaginary Y) of: 

Xe(R) + ^Gz(R,R) = i?. 



(24) 



(25) 



(26) 



The argument in the exponential is 2tt times the number of bosons enclosed by the imaginary momentum trajectory. 
To compare, in the field theory for the condensate we find the tunneling rate from: 



EwKB oc exp 



'Poe 



Wh -TT- drdr 

OT 



(27) 



where '0b(r, r) is the bounce solution of Eqs. k2 and E3. 

Whereas the semiclassical recipe for the tunneling of ij} is well defined, in the collective coordinate treatment it is 
not clear between which pair of trajectories the tunneling will take place, or equivalcntly, which value of E is to be 
used in Eq. |25|. This ambiguity stems from the fact that the collective coordinate approach does not really describe 
the vortex nucleation process, instead it describes the tunneling of a preexisting vortex bound to the disk in a closed 
trajectory towards an open trajectory further inside the flow. 

One can imagine doing the same kind of approximation used to obtain the formula for the critical velocity and 
simply admit that the tunneling will occur between the closed trajectory that reaches a distance 5 from the disk and 
its open counterpart. This would correspond to the vortex being produced at the furthest possible distance inside the 
compressibility dominated region around the disk. This fixes the value of E to be used in Eq. Pq and is consistent 
with the formula for the critical velocity, Eq. H^, since we then obtain Ewkb ^ 1 as Woo — > v^ (when the saddle point 
is exactly at a distance 5 from the disk). 

In Fig. ra we show the dependence of the tunneling rate with the applied velocity, Voo , calculated using Eq. BB for 
radii a = 1, 2 and 3^ and S — 0.5 and 1.5^. One observes, for small values of Woo/^cj that 



exp [-C [vc/vooT] 



a = 2.8 ±0.1, 



(28) 



is a good fit to all curves. 



VI. CONCLUSION 

We have analyzed the process of vortex depinning from a two dimensional obstacle due to an applied superflow. We 
discussed the process of dynamical depinning, above the critical velocity, and the process of depinning via MQT, at 
zero temperature and below the critical velocity, with a collective coordinate approach. The problems inherent to this 
approach were circumvented with the introduction of a length 5 {'^ ^) proportional to the size of the compressibility 
dominated region around the defect. A more thorough calculation of the critical velocity was made using the full 
NLSE and compared with the estimate based on the collective coordinate. The process of dynamical nucleation was 
demonstrated with a numerical simulation of the time dependent NLSE. 

It is left for a future work to contrast the the semiclassical tunneling rate computed from Eq. E5^ with the one 
obtained by solving the imaginary time NLSE. 

Experiments on critical velocities of two dimensional neutral superflows past obstacles do not exist in the literature, 
nevertheless one can expect the calculations presented here to be applicable to layered superconductors in the super- 
clean limit, when the pancake vortices are expected to follow the dynamics discussed in section III, see Ref. [p2i. When 
such systems are irradiated by heavy ions defects similar to the ones discussed here can be formed. Measurements of 
the critical current dependence with the temperature and applied magnetic field already exists, see for instance Ref. 
p3[ , but the dependence with the defect radius is still lacking, when obtained a comparison could be made with the 
results presented here. 

APPENDIX: KELVIN'S THEOREM 

Here we derive Kelvin's theorem, Eq. ^ from Eq. Q We consider a flow field v(r,t) and compute the rate of 
change of the circulation in a path C{t) being dragged by the flow: 



* v(r,<)-dL 
Jc(t) 



dt .ic( 



We must evaluate. 



1 



* v{r,t + dt)-d\~ f v(r,<)-dl 

Jc(t+dt) Jc(t) 



(Al) 



(A2) 



iC(t+dt) -ICit) 

The mapping between points in the C{t) loop and points in the C{t + dt) loop is 

r ^ r + v(r, t) dt. (A3) 

Therefore an element d\ in C{t) is mapped to 

r + dl ^ r + dl + v(r + dl, t) dt, (A4) 

d\^d\+ {d\ ■ V)v(r, t) dt. (A5) 



The integral in the C{t + dt) loop can then be written as 

* v[r + v(r, t) dt, t + dt] ■ {dl + {d\ ■ V)v(r, t) dt}, 
Jc(t) 



lc(t) 
or, to lowest order in dt. 



(A6) 



/ {v + ^t^r dt + {dt v ■ V)v} ■dl+ (f v • (dl • V)v dt. (A7) 

Jc(t) Jc(t) 



'C(t) JC(t 

Finally we obtain for the derivative 



d_ 
dt 



(f V • dl = / 9tv dl + / (v V)v • dl + / ivw^ • dl. (A8) 

Jc(t) Jc(t) Jc(t) Jc(t) 2 



The first two integrals add up to zero by virtue of Eq. H, whereas the third integral is identically zero as long as 
the path does not cross any singularity of w^ . 



The fact that the circulation is a constant along any path flowing with the fluid imply that the singular vortices 
must follow the flow velocity at their center. That is, Eq. nffl must hold. 
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FIG. 1. Isolines oi E = Xc + \Gi, see Eqs. y_l| and |l3|. These lines also represent the trajectories of a positive vortex in the 
presence of an applied flow with Ucx) = 0.4c along the y direction. The disk radius is a = 1^. The solid lines are below the 
saddle point, the dashed lines are above. 

FIG. 2. Critical velocity, in units of the speed of sound c, as a function of the disk radius a. The solid line was obtained from 
Eq. hTlwith 5 = 0.5^ and the dashed line from Eq. M. The squares are lower limits for Vc obtained from the time independent 
NLSE, Eq. Ell, with a grid spacing of 0.5^. The open circles used a grid spacing of 0.25^ (the two values coincide for a — 0.5^). 

FIG. 3. Countour lines of the condensate amplitude obtained from a numerical integration of the time dependent NLSE, Eq. 
[U The amplitude is zero inside the disk of radius 2.5^ and increases to the equilibrium value of y/po away from it. The initial 
state had a smooth circulation around the disk in addition to the applied flow in the y (vertical) direction of v^o = 0.5c. The 
times are 14, 23, 29 and 34^/c respectively for figures a), b), c) and d). In a) the vortex is beginning to form, in b) it is fully 
formed and being dragged by the flow, in c) one observes a vortex/ antivortex pair (the antivortex to the left and the vortex to 
the right) beginning to form and fully formed in d). 

FIG. 4. Close up view of the saddle point of Fig. hi Here the solid lines are the solutions of Eq. Eq with real Y and the 
dashed lines are solutions with imaginary Y. These can also be viewed as the real and imaginary time trajectories of a positive 
vortex around the disk of radius a = 1^ when an applied flow of Vaa ~ 0.4c along the y direction is present. The area enclosed 
by the dashed paths is proportional to the logarithm of the tunneling rate. 



FIG. 5. Logarithm of the semiclassical tunnehng rate as a function of the apphed flow velocity at infinity, calculated with 
Eq. ESl The value of E used in that equation was chosen to describe a vortex tunneling from a distance 5 from the disk towards 
an open trajectory, of same energy, further from it. Vc was calculated using this value of 5 in Eq. nTl The squares, circles and 
triangles correspond to disk radius a — 1, 2 and 3^ respectively. The solid symbols use S = 0.5^ and the open symbols use 
5 = 1.5^. 
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